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ABSTRACT 



An "RF sensor" network can monitor RSS values on 
links in the network and perform device-free localiza- 
tion, i.e., locating a person or object moving in the area 
in which the network is deployed. This paper provides 
a statistical model for the RSS variance as a function of 
the person's position w.r.t. the transmitter (TX) and 
receiver (RX). We show that the ensemble mean of the 
RSS variance has an approximately linear relationship 
with the expected total affected power (ETAP). We then 
use analysis to derive approximate expressions for the 
ETAP as a function of the person's position, for both 
scattering and reflection. Counterintuitively, we show 
that reflection, not scattering, causes the RSS variance 
contours to be shaped like Cassini ovals. Experimental 
tests reported here and in past literature are shown to 
validate the analysis. 

Categories and Subject Descriptors 

C.2.1 [Computer Systems Organization]: Computer- 
Communication Networks — Distributed Networks 

General Terms 

Design, Theory 

Keywords 

Fading, Device-Free Localization 

1. INTRODUCTION 

In this paper, wc develop methods for the analysis of 
wireless sensor networks which perform device-free lo- 
calization (DFL) of people and objects. DFL networks 
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locate people or objects without requiring them to co- 
operate in the system by carrying a radio device or a 
transmitting tag HOI LUl [HI [221 EOl [23 [23 ES] • Instead, 
DFL networks measure changes in the received signal 
strength (RSS) between many pairs of wireless sensors 
in the network to perform estimation of the locations 
of moving people and objects in the deployment area. 
This is advantageous in applications like emergency re- 
sponse, or secure buildings, in which one cannot expect 
that people who must be located will be wearing a trans- 
mitter tag. 

The DFL networks we study in this paper have five 
distinguising characteristics. First, they are device-free, 
in that the person (in this paper, wc use "person" gener- 
ally to refer to any large moving object) being tracked 
is not assumed to be wearing or attached to a trans- 
mitting or receiving radio device. Second, sensor nodes 
are radio frequency (RF) sensors. RF has particular 
advantages of being able to penetrate (non-metal) walls 
and smoke. Thirdly, they are sensor networks, in par- 
ticular, peer-to-peer networks of RF sensors, and thus 
measurements can be made between many pairs of RF 
sensors. Fourthly, they are model-based, i.e., sensor lo- 
cations are known, but there is no training nor a priori 
environment knowledge. Finally, they are narrowband. 
It is also possible to use ultra wide-band signalling to 
perform DFL [21 [2lJ [20] . This paper focuses on less ex- 
pensive narrowband sensors, which cannot distinguish 
multipath which arrive with different time delays. 

Variance of RSS due to human or object motion is due 
to changes in individual multipath which contribute to 
the received signal. The complex baseband voltage for 
a narrowband signal measured at a receiver is [4j, 
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where Vi = | 1 e- 7 * >i is the complex amplitude of the 
ith multipath component, and N is the total number of 
components. The received power is equal to \V\ 2 , the 
squared magnitude of the complex baseband voltage. 
What we call the "received signal strength" (RSS) is 



actually the received power in decibel terms, that is, 
R dB =201og 10 ||y|| 

Fading is the effect of the changes in phases {<j>i\i. De- 
pending on these phases, the sum in ([JJ may be destruc- 
tive (opposite phases) or constructive (similar phases). 
When objects and people in the environment move, they 
affect a subset of the multipath components depending 
on their current position. The new person or object 
may also scatter radio power and thus induce a new 
multipath component in the channel. In the bistatic 
radar literature, this new scatter-path is the only chan- 
nel change considered [17j . In this paper, we provide 
new models which allow us to quantify the total fading 
caused by changes in existing multipath components, as 
a function of the position of the new person. 

Fading due to human motion has been quantified in 
past studies in order to aid in the design of communi- 
cations systems which operate in such fading environ- 
ments pp. Researchers have observed that the motion 
of people near either a transmitter or receiver impacts 
measured fading statistics, and when more people are 
moving in proximity of a transmitter or receiver, fading 
increases [7j. 

Of particular interest for DFL, experimental results 
have been presented which relate a person's coordinate 
(with respect to TX and RX coordinates) to the vari- 
ance measured on a link. In [28], a contour plot, repro- 
duced in Figure [IJ shows that the "variation" (sum of 
absolute values of the differences of the time series) is 
highest when a person is located close to either node, 
is relatively high when a person is in the line segment 
between the nodes, and decays away from this line and 
the nodes. In contrast, experimental results in [30] show 
a different result, that the "dynamic" (average absolute 
value of the difference from the static mean) is highest 
in an oval centered at the midpoint of the line segment 
between the two nodes. We show in this paper how 
both results can be explained using analysis. We also 
present in this paper results of a new experiment con- 
ducted in our University bookstore which also validates 
the analysis. 

The goal of this paper is to provide a statistical model 
which relates (1) the variance of channel fading mea- 
sured on a link to (2) the position of the moving per- 
son in the environment w.r.t. the TX and RX locations. 
We desire a model which is based in the physics of 
radio wave propagation and does not require specific 
knowledge of the environment. For these purposes we 
apply traditional statistical spatial propagation models 
for scattering and reflection, which were originally de- 
veloped to study multi-antenna or directional antenna 
systems [TTI [To] . 

We obtain the statistical model in three steps. First, 




Figure 1: From [28j, fading signal "variation" as 
a function of a person's position. 

in Section [3J we show that the ensemble mean (across 
many realizations of a link) of RSS variance is approx- 
imately linearly related to the expected total affected 
power (ETAP) in dB. Secondly, in Section^ we present 
statistical multipath models for scattering and reflec- 
tion. Finally, in Section we quantify the ETAP a 
function of the person's position (with respect to the 
TX and RX positions). 

A counter-intuitive result of the analysis is that RSS 
variance has contours similar to Cassini ovals, as would 
be seen in a bistatic scattering power equation; however, 
this dependence occurs when the sole multipath propa- 
gation mechanism is reflection. We verify the new model 
using measurements in Section [7]- those reported in the 
literature [28[ 130] , and an extensive new measurement 
campaign we conduct and present in this paper. 

2. RELATED WORK 

In the Introduction, we reviewed propagation research 
which reported the effects of human motion on fading 
statistics. DFL takes advantage of advances in wireless 
sensor networks to develop networks of "RF sensors" for 
which the radio is the sensor. Woyach, Puccinelli, and 
Haenggi [27J were the first to present RSS-based DFL, 
to our knowledge, and termed these applications as "sen- 
sorless sensing". As one example, they showed how RSS 
measurements from sensors deployed in a hallway could 
be used to track someone walking through, and to in- 
dicate whether they were walking or running. Zhang 
et. al. [301 [SI] positioned sensors on a ceiling and used 
the change in RSS on links to estimate the location of 
people moving below. 

Youssef, Mah, and Agrawala [29] coined the term "device- 
free passive (DfP) localization" to distinguish this sys- 
tem, in which people are located regardless of whether 
they carry a radio, from active RFID systems which lo- 
cate a TX or allow a RX to locate itself. Two detection 
methods and a radio map-based tracking method are 
presented in [25]. Other work has also used training 
measurements to perform localization and tracking [231 
1221 114] . Because training measurements require signif- 



icant time, automatic generation of radio maps is "one 
of the main technical challenges behind the implemen- 
tation of DfP" and it is critical to know "the effect of 
these entities [in the environment] on the received signal 
at different points in the environment" [29] . This paper 
provides statistical models and tools for this purpose. 

Our past research has also used changes in the statis- 
tics (mean and variance) of RSS in order to estimate 
an image of presence of people in the environment [24j 
125] [26] . In particular, we have shown the ability of an 
RF sensor network to track a person moving inside of a 
building through external walls [26] . The models devel- 
oped in this paper are useful for development of estima- 
tion bounds and improved estimation algorithms. 

Other related efforts have explored details regarding 
the effect of human movement on measured RSS. An 
experimental test in [10] uses RSS on a link to reliably 
detect the presence of a single person walking from one 
node to the other, in a variety of different indoor and 
outdoor environments. Experiments in [15j go further 
by providing linear relationships between the RSS mean 
and variance and the number of people walking or sit- 
ting between two nodes, for use in crowd density es- 
timation. Very relevant to this paper, measurements 
reported in [28] show the "variation" of RSS caused by 
human movement as a function of position across the 
2D plane containing the two nodes. 

Research has also investigated polarization [18] as an 
alternative to RSS. When two antennas can be arranged 
to have orthogonal polarization, it adds significantly to 
the ability of a detector to measure a change in the 
environment [18j . 

Further, ultra wide-band (UWB) measurements have 
long been suggested for through- wall surveillance (TWS). 
Multi-static UWB approaches use networks of UWB 
transceivers to locate the position of a new person in the 
environment within and around the network and thus 
are closely related to this research [21 [17j [20] [21] . From 
the radar perspective, a new person in the environment 
creates an additional scattering multipath component. 
In multipath-rich environments, such as indoor environ- 
ments, the changes to existing multipath are at least as 
important as any additional scatter-path caused by the 
presence of a new person. This paper focuses on the 
effects of a new person on the existing multipath in the 
channel. 

3. RELATION OF RSS VARIANCE AND ETAP 

In this section we show that the variance of RSS has 
an approximately linear relationship with the expected 
total affected power (ETAP), that is, the expected value 
of the powers in affected multipath components. In 
other words, a person's new position causes random 
changes in some of the multipath components, and the 



sum of the power in these affected components has a 
linear relationship with the variance of fading. We start 
by showing how random changes in affected multipath, 
and no changes in other multipath, lead to a Ricean 
distribution of fading. We present a linear relationship 
between the variance of RSS (in dB) and the K fac- 
tor (in dB) of the Ricean distribution. Then we show 
that the ensemble mean of the K factor (in dB) has an 
approximately linear relationship with the ETAP. 

To define ETAP, consider the changes in the complex 
baseband voltage in {1} that occur when a moving per- 
son now exists at position x D . We assume that some of 
the multipath components i are affected in phase and 
amplitude. Other multipath components, because their 
path does not come close to the person, are unaffected 
by the change. Define sets T and T as the indices of 
multipath which are unaffected, and affected by the new 
person, respectively. ETAP is defined as 



ETAP = E 



E 

LiST' 
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(2) 



where Vi is the complex baseband voltage of multipath 
i. ETAP is a function of that person's position x Q , as 
will be discussed extensively in Section [5] but first, we 
show its relationship with RSS variance. 

3.1 Ricean Distribution of Fading 

As discussed above, we can consider multipath to be 
either affected or unaffected by the moving person at 
position x D . Consider partitioning the sum in © into 
two sums, 



(3) 



The first sum over i G T is unaffected by the place- 
ment of the new person, and therefore is constant (i.e., 
non-random) for a particular link. We denote this first 
sum as V — X^ier^- The second sum over i (£ T is 
random because the multipath components change ran- 
domly because of the motion of the person. We assume 
that the phases are randomly altered (uniform between 
and 7r) and multipath components i (£ T are multiple. 
Applying the central limit theorem, the real and imagi- 
nary parts of the sum of the affected multipath complex 
baseband voltages are seen to be independent and iden- 
tically distributed (i.i.d.) zero-mean Gaussian random 
variables V ns j and V ns ,Q [!]■ So, 



V = V + V ns j+jV nSiQ 



(4) 



The variance of V ns = V ns j +jV ns ^Q is equal to the sum 
of the magnitude squared of the voltages in its sum, that 

is, Zitr\Vi\ 2 m : 

The envelope ||V|| is Ricean because it is the magni- 
tude of a complex Gaussian random variable with non- 



zero mean V. If there were no unaffected multipath, i.e., 
|T| =0, then the envelope would reduce to a Rayleigh 
random variable. The if-factor of this Ricean distribu- 
tion is given by 



K = 



E 



ieT' 



(5) 



3.2 Relation of RSS Variance and K Factor 

In past research, we have shown that the variance of 
RdB i when the envelope is Ricean, is purely a function of 
the Ricean K-factor of the envelope [26]. In fact, there is 
an approximately linear relationship between the vari- 
ance of RdB, Var[i?ds], and the K-factor in decibels, 
KdB = 101og 10 if, for a range of —2 < KdB < WdB 
|26j . In this linear region, Var [RdB] varies from 27 dB 2 
to 3 dB 2 . This linear region is when the variance of 
RdB is most important. Stationary links will measure 
Var [RdB] < 3 dB 2 due to noise. In our experience, 
very few links with TX and RX stationary will measure 
Var [RdB] > 27 dB 2 . Thus, for our purposes, it is accu- 
rate to describe the relationship between Var [RdB] and 
KdB as linear. Furthermore, it is decreasing - as KdB 
increases, it means less severe fading is occurring, and 
thus a lower Var [RdB]- To be explicit, 



Var [R, 



dB 



a o ~~ a\K 



dB-. 



(0) 



for some positive constants ao, a\. 
In dB terms, ([5]) becomes 

K dB = 101og 10 [V] 2 - 101og 10 (]T |V5I 3 ) 

Effectively, KdB is the power in the sum of unaffected 
multipath minus the power in the affected multipath. 
What is the expected value of KdB, taken across an 
ensemble of environments with the same geometry of 
TX, RX, and person position? 



E[K dB ] =£[l01og 10 |^| 2 ]-£ 



10 log- 



in X 



Vi\ 



(7) 



In relative terms, the power in the unaffected multi- 
path are generally much stronger than the power in the 
affected multipath. With multipath dispersing in all di- 
rections around an environment, a person typically af- 
fects only some small fraction of the total power (unless 
the person is right on top of the TX or RX) . If we rewrite 
the left term, E [10 log 10 [V] 2 ^ , as the expected value of 
the 101og 10 of (total power - affected multipath power), 
we can see that it is largely constant when the affected 
multipath power is much less than the total power. We 
approximate ([7|) as, 



E [K d B] ~ c up - E 



lOlogio^l^l 



(8) 



for a constant c up . When the affected power is one half 
of the total power, this constant approximation is in 
error by 3 dB. When the affected power is less, the error 
is respectively less. 

Combining ((5|) and ©, we have that 



E [Var [R dB ]] 



a,2 + a\E 



101og 10 w\ 



iET> 



for ai = ao — a\c up . Note again that the variance is 
taken in one particular environment for one particular 
link, while the expected value is across the ensemble of 
environments and links with the same geometry 

Finally, we approximate E [log Y] ss log E [Y] , for Y = 
J2ieT' l^il 2 - We show in Appendix [El that this is a good 
approximation when |T'| is more than one or two, i.e., 
multiple multipath are affected by the person's position. 
With this approximation, 



E [Var [R dB ]] w a 2 + ffll 10 log 10 E 



En 



(9) 



In summary, the ensemble mean of RSS variance has 
an affinc relationship with the expected value of the sum 
of the powers in affected multipath components, i.e., 
ETAP. Sections H and [5] analyze ETAP, in particular, 
how it is a function of the propagation mechanisms and 
TX, RX, and person's positions. 

4. STATISTICAL MULTIPATH MODELS 

To quantify fading as a function of an person's posi- 
tion in space, we need a statistical channel model which 
describes the spatial extent of each multipath in the 
channel, what we term a spatial multipath model. Sev- 
eral spatial multipath models have been developed in 
past research for the study of directional antennas in 
cellular systems [TTJ [TB], for the development of fad- 
ing models [3], and for multiple-input multiple-output 
(MIMO) systems [5], among other purposes [5]. We fo- 
cus on models which explain the geometric path of each 
multipath component, rather than those based purely 
on measurements which model only angle-of-arrival (AOA) 
or time-of-arrival (TOA) at the receiver. We also need 
models applicable to the case where antenna heights are 
both relatively close to the ground (as opposed to being 
at cellular base station heights). 

4.1 Related Work 

The geometrically-based single-bounce model (GBSBM) 
jll| is a spatial multipath model used in cases when 
the antenna is at the same height as the scatterers. 
In this model, each (non-LOS) multipath experiences 
a single reflection. Any multipath component which ar- 
rives with excess delay less than a threshold is known to 
have changed direction at some point within an ellipse. 



These reflection points are assumed to be uniformly dis- 
tributed in the environment. In the received power 
at location x r of the multipath component originating 
at x t and reflected at x is given by, 



P,(x) = 



! x « - x l! 



(10) 



where c r is a constant, and n p is the path loss exponent. 

The model of [IS] is also a geometrical model which 
assumes each path is scattered once by a scatterer on 
a plane. In |16j . however, the plane does not contain 
the TX and RX antennas, which may exist at arbitrary 
heights }ibs and huts for the base station and mobile 
station, respectively. The model of [T^] is also single- 
bounce, but the propagation mechanism is scattering, 
so that the received power of a scattered multipath is 
given by, 



^a(x) 



x 



(11) 



where c s is a constant. 

The locus of points {x : P s (x) = 7} for some constant 
7 is called a Cassini oval with foci at Xt and x r |17j . For 
different values of 7, Cassini ovals are plotted in Figure 
[2J In general, P s (x) is low when the scatterer at x is far 
from both the TX and RX, and is high when x is close 
to one of them or the line in between them. 
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Figure 2: A contour plot of Cassini ovals for the 
case when TX and RX locations are shown as ■, 

and c, = 100. 



4.2 Our Model 

For the statistical multipath model in this paper, we 
consider the union of the models of [TT] and [16] . To be 
specific, we make four simplifying assumptions about 
the radio channels in the environment of interest: 

1. Omni-dircctional antennas are used at each node. 



2. Mechanisms for propagation include scattering and 
reflection. Diffraction is not considered in this pa- 
per. 

3. The scattcrcrs (and reflectors) are located in a plane 
parallel to the ground, the "scatterer plane". 

4. Single-bounce propagation. Each multipath has a 
single change in direction on its path from TX to 
RX. 

As they were in [11] and |16| . these are significant sim- 
plifying assumptions. The single-bounce model is, as 
Chang and Sahai concisely described, "a bit dubious" 
[2], but it can help to provide a first-order model for 
multipath- induced characteristics. More general anal- 
ysis which assumes multiple reflections and directional 
antennas can be an extension of techniques described 
in this paper. Diffraction loss has a complicated rela- 
tionship with object geometry and angles of arrival and 
departure from the object, which makes general analysis 
difficult. 

When we consider multipath which have experienced 
a single scatter at location x, we describe their received 
power as P s (x) as given in ([lip . When we consider 
multipath which have experienced a single reflection at 
location x, we describe their received power as P r (x) as 
given in (j!0[) . To keep notation consistent, we describe 
the physical item which exists at location x 
terer" regardless of whether the item reflects or scatter- 
ers the radio wave. Each different scatterer might have 
a different constant c s or c r , however, so for simplicity, 
we assume that both arc constant across the ensemble of 
scatterers. We have assumed isotropy of antennas and 
radar cross-section (RCS) and thus c s and c r arc not a 
function of position. 

We also note that the power of a scatter-path from the 
new person can also be described using our model. In 
past research [T7J [301 125 it is assumed that a new per- 
son in the environment creates an additional multipath 
component based solely on scattering from the person. 
The received power of this new scatter-path from the 
new person at location x D can be written as P s (x D ). 

Similar to [TTJ [TB] we assume scatterer locations are 
Poisson distributed across the scatterer plane. W.l.o.g, 
the scatterer plane z-coordinate is taken to be zero. We 
denote the average scatterer density as ij/m 2 . The plane 
containing all RF sensors is also parallel to the ground, 
i.e., the TX and RX locations are always at the same 
height, which we denote Az. 

We also assume that the person can be approximated 
by a vertical cylinder as shown in Figure [3] In the scat- 
terer plane, it is a circle centered at x Q with diameter 
D. The top of the cylinder is higher than the highest, 
and the bottom of the cylinder is lower than the lowest, 
of Az and zero. This is only a crude approximation of 



a person's shape, but is simple to analyze because the 
cylinder casts an infinite "shadow" onto the scatterer 
plane. Given a scatterer at location x, if either the line 
segment connecting the TX to x, or the line segment 
from x to the RX, intersects with this vertical cylinder, 
then the multipath is "affected". 
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Figure 3: TX, RX, plane containing scatterers, 
and the new person. 



5. ANALYSIS 

Given the model assumptions described in Scction|421 
we analyze two types of ETAP, first assuming the prop- 
agation mechanism is scattering, and then assuming it 
is reflection. Certainly, any real-world link will have 
multipath caused by both reflection and scattering. In 
addition, there may be a line-of-sight path which is not 
reflected or scattered. The real-world ETAP will be 
some linear combination of the ETAP due to scattering 
and the ETAP due to reflection. 

Key to this analysis is to see that single-bounce mul- 
tipath components affected by a person have scatterers 
in the shadow of that person, with respect to cither 
the TX or RX. By "shadow" with respect to the TX or 
RX, we mean that the line segment from either the TX 
or RX to any point in the shadow crosses through the 
person. Figure |4] shows a multipath component with a 
scatterer in the shadow of the person with respect to 
the RX. The analysis of affected multipath transforms 
into an analysis of the positions of scatterers within an 
area A = A t U A r , where we denote the areas shadowed 
by the person with respect to the TX or RX as A t and 
A r , respectively. 

The area A t is essentially an infinite height isosceles 
trapezoid; equivalently, an infinite height isosceles trian- 
gle with the area between the vertex x t and the person 
at x D removed. Figure [4] shows the median line, i.e., 
the line from x t which passes through x D . We make two 
simplifying assumptions: 

• Approximation 1 : The width of the shadowed area 
at x Q is equal to D, the diameter of the person. 




Figure 4: The TX x t , RX x r , and the "affected" 
area (shaded) in which a scatter could be lo- 
cated in order to have a single-bounce path that 
crosses the person located at x D . 

• Approximation 2: The value of a function /(x) for 
x in the shadowed area can be approximated by 
/(x), where x is the projection of x onto the me- 
dian line, i.e., the point on the median line closest 
to x. In other words, the isosceles trapezoid of the 
shadowed area is narrow enough that the value of 
/(x) can be assumed to be nearly constant along 
the perpendicular direction from the median line. 

Both approximations are good when D -C ||x — x t ||. 
These approximations are also used with A r , the shad- 
owed area by the person with respect to the receiver, by 
replacing x t with x r in the above description. 

In a homogeneous Poisson spatial scatterer process, 
scatterer locations arc randomly distributed across the 
scatterer plane with density r\ scatterers per unit area. 
In each ETAP (scattering or reflection) , we want to find 
the expected value of a function of the scatterer loca- 
tions. The ETAP is denoted Q and is written as, 



Q = E 



/(Xi 



where the sum is over every scatterer location, and where 
/(xj) is a scalar function of scatterer position Xj. In 
general, we write this expected value as, 



Q = 



[ [ 77/(x)dx. 



(12) 



At this point, we make another approximation: 



• Approximation 3: We approximate the area inte- 
gral in (|12[) by the sum of two shadowed area in- 
tegrals, one over A r and the other over A t . With 
this approximation, Q = Q t + Q r , where 

Qt= 1 1 rjf(x)dx, Q r = t]f(x)dx 

(13) 

Essentially, we assume that the two parts of the 
shadowed area, A r and A t , are nearly disjoint. 

This approximation is good except for when x D is on the 
"far" parts of the line containing x f and x r , i.e., on the 



straight line containing x t and x r but not in between 
Xj and x r . On these far parts of the line, the shadows 
of the person w.r.t. the TX and RX overlap, and we do 
not consider the resulting expression to be accurate. 

Because of approximations 1 and 2, the area integrals 
in (fl~3|) can be simplified to be written as line integrals 
as follows. First, consider the area shadowed w.r.t. the 
TX, A t . Because of approximation 2, the value of /(x) 
is constant along the line segment perpendicular to the 
median line within the shadowed area. To simplify, we 
write x as the projection of x onto the median line. The 
vector x can be written as, 



for some a > 0. The fraction (x D — xt)/||x — xt|| is 
the unit vector parallel to the median line. The width 
of the shadowed area at x is the length of the line seg- 
ment perpendicular to the median line contained with 
At . This width is calculated using similar triangles. The 
isosceles triangle that starts at x t and ends at the per- 
son location x Q has height ||x< — x Q || and has base length 
approximately equal to D (by approximation 1). The 
full isosceles triangle has height ||x t — x D | + a, and we 
denote its base length as b t . So, 



h = D 



a + ||x t - x G | 



||x t - x || 
Thus (fT3|) simplifies to become 

Q = Qt + Qr, where (14) 

[°° a + ||xj - x Q || / Xg-Xj \ 
Qt = / r)D — n n — f\x + a T [7 ) da 

Qr = / V D — n n — / I x o + a n 1, da 

Ja=0 \\ x r x o|| V ll x o — x r||/ 

5.1 ETAP: Scattering 

First, we consider the expected total affected power 
(ETAP) assuming that all multipath are due to scatter- 
ing. We use the scattering power formula of (fTTj) as the 
function /(x) in (fH)) . 

We note that there is a singularity in (fTTj) whenever 
the scatter position x is equal to the TX or RX position. 
Infinite received power is not physically possible - it is 
an artifact of a far-field approximation formula evalu- 
ated at near-field - thus it is appropriate to ensure that 
P s (x) is finite for all positions x. The singularity can be 
remedied by adding a constant into the denominator, as 
studied in [8]. Instead, we have specified that the scat- 
terer plane does not contain the TX or RX locations, 
thus avoiding the singularity. 

We apply /(x) = P s (x) to find Q = Q t + Q r - For Q t , 



the integral in (|14| simplifies to, 

r jwnx t -x | 

' a=0 (Hxt-xoll+a) (x r -x G )- 



Qt = 



l|Xo-X t || 



; da 



(15) 

The expression for Q r is obtained by switching x t and 
x r in (Ti"5|) . We show in the Appendix that Q simplifies 
to 



Dc s r] ( (n - 0)(1 + cos( 



d 2 

u rt { u x 



sin 9 



— \o H x * Z x °l 
d Xn ||x r -x | 



where d rt = ||x r — x t ||, and 

d t = (||x r - x || _1 + ||xj - x || _ 

rf Xo = (ll x '' - Xoll -1 - ll x < - x o|T' 

'(x r - x ) T (x - x t ) 



(16) 



cos 



x - x t 



(17) 



Note that d+ is the combination in parallel of the two 
distances from the TX and RX to the person, and d~ 
is the parallel difference of the two distances. The angle 
8 is the angle between the line from x Q to x r and the 
line from x t to x . Note that because of the cos _1 (-), 
< 9 < it. 

For the case of x t = [—1, 0, 0.1] T and x r = [1, 0, 0.1] T , 
the ETAP as a function of x G is shown in Figure [5ja). 
Note that since the plot is shown in dB relative to the 
maximum value, the constants D 1 c s , and rj do not 
change the plot. 

5.2 ETAP: Reflection 

In this section we calculate the ETAP assuming that 
reflection is the only propagation mechanism. In this 
case, the function /(x) is set equal to P r (x) as given in 
(|T0|) . We still have that Q = Q t + Q r , and the integral 
expression for Q t is 



Qt = 



" ULr ||x t -Xo|| 



Q=0 



a 



a-; 



-x t || 



nr(18) 



Unfortunately, this Q t is not generally tractable for ar- 
bitrary real-valued path loss exponents n p . Instead, 
we perform numerical integration to find Q t and Q r . 
Compared to the exact equation for reflection ETAP, 
the expression in (|18[) does simplify calculation by re- 
quiring only one integral. For the case oi n p = 3, 
x t = [-1,0,0.1] T and x r = [1,0,0.1] T , the ETAP as 
a function of x Q is shown in Figure [51(b) . Again, since 
the plot is in dB relative to the maximum, the constants 
D, c r , and r\ do not change the plot. 

5.3 Parameter Sensitivity 

In this section, we study the effects on ETAP of changes 
in path loss exponent, n p , and the relative height of the 
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Figure 5: ETAP in dB relative to the maximum, 
with xt and x r as shown (■) and Az = 0.1m, for 
(a) scattering, and (b) reflection (with n p = 3). 



by providing greater location precision when ETAP (or 
variance of RSS) is measured to be high. 



-10 



-15 



-20 



g -25 



-30 



-35 



-40 



" 4 - 5 2 



























-■-2.0 






3.0 






- - 4.0 
5.0 











X Coordinate 



Figure 6: Reflection ETAP (dB max) 

x Q with y-coordinate equal to 0.1, for n 

{2.0,3.0,4.0,5.0}, and x t = [-1,0,0.1] T and 
[1,0,0.1] T . 



for 

p £ 
x r = 



Next, we consider the effect of the separation between 
the TX/RX plane and scatterer plane, Az. For the case 
of reflection, the effects are shown in Figure [7]Jb) . For 
low values of Az, when the scattering plane contains or 
nearly contains the TX and RX, there are noticeable 
peaks in ETAP near the TX and RX. In contrast, when 
the TX and RX are in a plane high above (or low be- 
low) the scatterer plane, the shape of the ETAP surface 
actually turns around such that its maximum is halfway 
between the TX and RX. 

For the case of scattering, the effects are shown in 
Figure [7fa). In the scattering case, the ETAP is always 
highest halfway between the TX and RX. But, with 
increasing Az the shape of the ETAP becomes rounder, 
and the slope outside of the area between the TX and 
RX becomes shallower. 



plane containing the TX and RX, Az. 

For the reflection ETAP, we vary the path loss expo- 
nent n p . We show the results by plotting in Figure [6] 
a horizontal cut of the value of ETAP (along the line 
with y-coordinate equal to 0.1). Note that all plots are 
in dB relative to the maximum of ETAP across all space 
and for all tested values of n p . We observe that increas- 
ing path loss exponent generally decreases the ETAP. 
Further, away from the TX and the RX (below -1 and 
above +1), an increasing n p also increases the rate of 
decrease in ETAP. In effect, the with high ETAP 
is more isolated from the area with low ETAP when 
n p is higher. So although higher n p reduces the mag- 
nitude of the ETAP, it may benefit radio tomography 



6. DISCUSSION 

A few basic assumptions and approximations have al- 
lowed analysis of affected multipath, as a function of 
x D and the TX and RX positions. These quantities in- 
clude the total number of affected multipath and the 
ETAP when propagation is due to scattering or reflec- 
tion. Previous research has surmised that the new mul- 
tipath caused by scattering from the person was the sole 
change in the channel; and that this change introduced 
a new path with power given by the bistatic scattering 
equation in ()11|) . The contours of this power surface 
P s (x ) have the Cassini oval shape. 

The ETAP expressions in this section, in combination 
with the ETAP/RSS variance relationship presented in 
Section [3l show the "shape" of the typical RSS vari- 
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Figure 7: ETAP (dB max) vs. x Q with y- 
coordinate equal to 0.1, for several Az, when 

x t = [-l,0,Az] T and x r = [l,0,Az] T , for (a) re- 
flection and (b) scattering. 



ance as a function of person position x Q . By considering 
ETAP, in contrast to past research, we have shown that 
reflection, not scattering, leads to RSS variance surfaces 
with contours similar in shape to Cassini ovals. The re- 
flection ETAP has two peaks around the TX and RX 
positions, except when Az is large. In contrast, Figure 
03a), which shows the scattering ETAP, does not show 
peaks around the TX and RX positions for any Az. 
Instead, Figures [5{a) and EOa) show a single rounded 
rectangular peak when x D is between the TX and RX 
locations. 

In summary, when reflection is the dominant mech- 
anism, and the scattering plane is close to the plane 
containing the transceivers, the maximum effect is seen 
when the person is very close to either the RX or TX. 
When scattering is the dominant mechanism, or when 
the scattering plane is far from the plane containing the 
transceivers, the maximum effect is seen when the per- 
son is halfway between the RX and TX. 



7. EXPERIMENTAL RESULTS 

We refer in the Introduction to two studies which re- 
ported experiments which quantify the RSS variation as 
a function of the TX, RX, and person's location. The 
results of [55], shown in Figure[TJ show highest variation 
when the person is close to either node and high varia- 
tion when the person is close to the line connecting the 
nodes. In contrast, experimental results in [30] Figs. 4, 
5] show a different result, that the "dynamic" (average 
absolute value of the difference from the static mean) is 
highest in an oval centered at the midpoint of the line 
segment between the two nodes. 

The differences in the experimental setups may ex- 
plain the different results. In [28], RF sensors are lo- 
cated on tripods, at an apparent height of about 1.5 
meters, thus the plane containing the TX and RX is 
likely to also include many scatterers. Since Az is close 
to zero the analysis would tell us to expect the highest 
lines from Figures E^a) and (b) . Because the variance 
is highest closest to the TX and RX, the analysis then 
suggests that the main propagation mechanism is reflec- 
tion. 

In contrast, in [30] . all RF sensors are located on the 
ceiling (at 2.4 meters) in an empty room. Because the 
room is empty, the scattering plane (likely the floor) is 
separated by Az = 2.4 meters from the plane containing 
the TX and RX. Neighboring sensors are separated by 
2 meters. Since Az is greater than the distance between 
sensors, the analysis indicates to expect the lowest lines 
in Figures [7]J a) and (b). Both reflection and variance 
ETAP at high Az have a shape which is highest in the 
midpoint between the TX and RX, similar to the exper- 
imentally reported results. 

We perform an extensive experiment using the fol- 
lowing procedure to experimentally determine the RSS 
variance surface in another environment. Thirty-four 
nodes (Crossbow TelosB, operating the IEEE 802.15.4 
protocol at 2.4 GHz) are placed throughout a 16 by 36 
feet area within the University of Utah Bookstore, all 
at approximately 1.2 meter height. It is a typical store 
environment, with bookshelves, books, aisles, tables, 
chairs, product displays, and other obstructing objects 
at height similar to the node height. The location of 
each node is surveyed and recorded. The signal strength 
of each link is recorded over time as a human walks in 
a known track within the area. For each measurement 
of RSS, RdB, on each link, the coordinates of the TX, 
RX, and human are all known. We rotate, translate, 
and scale the coordinates such that x t = [1, 0] T and 
x r = [— 1,0] T , that is, we normalize the coordinates 
so that the TX and RX coordinates are the same as 
used in the examples in Figure [5] The person's posi- 
tion is transformed in the same manner to obtain x Q so 
that its relative position w.r.t. x t and x r is preserved. 
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Figure 8: Experimental RSS variance vs. x Q , nor- 
malized such that the RX is located at (-1, 0) and 
the TX is located at (1,0). 

In post-processing, we create a two dimensional grid of 
bins across the range of x Q , and each RSSI measurement 
RdB for each time sample is stored in a list for the bin 
that corresponds to the human position x D . This pro- 
cess is repeated for every link measurement made during 
the experiment, over 110,000 measurements. The vari- 
ance of RSS measurements for each bin is calculated 
and plotted as a surface. A contour plot of the variance 
surface is given in Figure [8j 

Because in our experiments the scattering plane is 
similar in height to the plane containing the RF sen- 
sors, we have Az close to zero. Similar to [28], we by 
comparing the results in Figure [5] to the ETAP analysis 
that the dominant propagation mechanism was reflec- 
tion. 

8. CONCLUSION 

This paper presents models for the variance of RSS 
caused by motion of a person (or object), as a function 
of the person's position, x D , and the TX and RX posi- 
tions, xt and x r . The analysis is statistical, rather than 
based on any specific knowledge of the environment, and 
considers single-bounce multipath caused by scattering 
and reflection. We expect to see surfaces with contours 
of Cassini ovals when considering scattering from the 
person. However, this new analysis shows two interest- 
ing things: 

• Even when multipath are all caused by reflection, 
variance surfaces vs. x„ can have contours similar 
to Cassini ovals. 

• When multipath are caused by scattering, or in 
high Az cases of reflection, variance surfaces are 



highest when the person is in the middle of the 
line between the TX and RX. 

We show using reported measurements, and a new ex- 
tensive measurement set, that the analysis shows results 
that are similar to those observed experimentally. 

Future analytical and experimental work is required 
to verify and extend the results presented in this pa- 
per. Analytical results should be extended to the case 
of multiple 'bounces', i.e., when multipath may expe- 
rience more than one reflection or scattering. Exper- 
imental work should measure the total affected multi- 
path power in a variety of controlled and uncontrolled 
environments to test the relationships derived in this 
paper. Finally, we will apply the developed statistical 
model to improve estimation bounds and imaging and 
tracking algorithms for use in DFL systems. 

APPENDIX 

A. PROOF OF SCATTERING ETAP 

In this section, we prove (|16p from (| 1 5[) . Expanding 
from ([T5j). 

da 



Qt = 



Dc s ri 

||x t - x || J a=0 (a + a)(a 2 — fa + b 2 ) 



where we define a 



and 



/ = 2||x r — x || cos6. By partial fraction expansion, the 
fraction in the integral can be written as, 

1 -|(2a-/) + i(2a+/)- 
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a 2 + af + b 2 
Note that 
a 2 + fa + b 2 = 
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cos 9 
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= || (x - x 4 ) + (x r - x )|| 2 = ||x r - x t | 
Using an integral table [5J 2.103(4)] we can obtain, 
„ Dc s r] 
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B 



2a + f 



log- 



arctan ■ 
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2a -f 
V^ 2 - I 2 



■\J 'a 2 — fa + b 2 

where arctan(x) is the arctangent of x. Note that [B|^L 
log 1 — log f . For the A term, we have that \J Ab 2 — f 2 = 
2||x r — x || Vl — cos 2 9 = 2||x r — x || sin (9, where the ab- 
solute value around the sine term is not necessary since 
< 9 < ir. Finally, 2a + f = 2(||x t - x || + ||x r - 
x || cos 9). Evaluating the limits on A, 
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And thus expressions for Q t , and similarly, for Q r , are: 
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(19) 



Plugging the final expressions for Q t and Q r into (fl4)) . 
we have the expression in (fT| 



B. APPROX. EXPECTED LOGARITHM 

In this section, we justify that E [10 log 10 ^1 ~ 10 log 10 E 
for Y = J2ieT' l^i| 2 - Converting to natural logarithms 
and representing \ Vi\ 2 = V?j + V?q, where Vij and V^q 
are the real and imaginary components of Vj, 



E[mog 1Q Y] = 



10 



log(10) 



log £ + 



(20) 



We can't evaluate (|2"D|) without making a distribu- 
tional assumption about Vij and V^q. It is typically 
assumed that {Vi}i are zero-mean complex Gaussian 
random variables [4]. In our case, we assume that V,j 
and Vi,Q are i.i.d. with variance a 2 . Clearly, not all mul- 
tipath components have the same variance, since mul- 
tipath caused by more distant scatterers tend to have 
lower power. However, the equal variance case is one 
that is analytically tractable and can be used to show 
the accuracy of the approximation. 

From (|20p , we renumber the multipath as 1 through 
2m, w.l.o.g, where m = \T'\, 



£[101og 10 Y] 



10 



log(10) 



-E 



log^l^ 

3=1 



(21) 



where Uj are i.i.d. Gaussian with zero mean and vari- 
ance a 2 . Moser [13] has shown that the expected value 
of the logarithm of this sum evaluates to 



£[101og 10 F] 



10 



log(10) 



'IE 



1m- 1 -, \ 



(22) 



where je is Euler's constant, 7# ~ 0.577. In this case, 
our approximation would be 

2m 

101og 10 £[F] = 101og 10 E£[|^| 2 ] = 101og 10 2ma 2 

3=1 

(23) 

We compare (j2"B"f and (|2^|) numerically. For one multi- 
path, the worst case, the approximation is in error by 



1.2 dB. With two multipath, the error is 0.6 dB, and 
with five or more multipath, the error is 0.2 dB or less. 

In summary, for the case of i.i.d. zero-mean com- 
plex Gaussian multipath, the error in approximating 
i£ [101og 10 y] as 10 log 10 E [Y] is small even for small 
numbers of multipath. As more multipath are present, 
the approximation becomes more exact. 
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